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Abstract 

An absolutely convergent double series representation for the density of the supremum of a- 
stable Levy process is given in [3, Theorem 2] for almost all irrational a. This result cannot be 
made stronger in the following sense: the series does not converge absolutely when a belongs to 
a certain subset of irrational numbers of Lebesgue measure zero (see [6, Theorem 2]). Our main 
result in this note shows that for every irrational a there is a way to rearrange the terms of the 
double series, so that it converges to the density of the supremum. We show how one can establish 
this stronger result by introducing a simple yet non-trivial modification in the original proof of [3, 
Theorem 2]. 
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1 Introduction 



Consider an a-stable Levy process X, defined by the characteristic exponent 

^{z) = -InE [e''^'] = 1^1" X (e™(5-p)i^^>o^ + e-™(^"'')l{,<o}) , zeR. (1) 

Here (a, p) belong to the set of admissible parameters 

A = {ae (0, 1), p e (0, 1)} U {a G (1, 2), p e [1 - a-\ a'^]}. 

It can be shown (see [6] and Section 4.6 in [8]) that with this parameterization we have p = P(Xi > 0) 
and that any stable process with a ^ 1 can be rescaled so that it's characteristic exponent is written in 
the form (1). 

Our main object of interest is the density p{x) of the supremum 5*1 := sup{Xu : < m < 1}. Let us 
summarize what is currently known about this function. 

(i) The simplest case is when X belongs to one of the Doney classes Ck,i (here k,l G Z), which are 
defined by the relationship p + k = I /a (see [2]). In this case an absolutely convergent series 
representation for p{x) is given in [5, Theorem 10]. From now on we assume that X does not 
belong to one of the Doney classes. 

(ii) When a is rational, [6, Theorem 3] gives an explicit formula for the Mellin transform of p{x), which 
is expressed in terms of elementary functions and the dilogarithm function Li2(^;). Unfortunately, 
there seems to be little hope of obtaining an explicit series representation for p{x): This problem is 
equivalent to evaluating the residues of the Mellin transform, which in this case has multiple poles 
and computing these residues in closed form seems to be impossible. 

(iii) When a is irrational, we define sequences {am,n}m>o,n>o and {bm,n}m>o,n>i as follows 

(-1)"^+" yr jap + j - 1)) A sin(7ra(p + j - 1)) 

T {1- p-n-^)T{ap + m + an) sin (^) "Li sin(7raj) ' 



r{l-p-n-^-^)T{ap + m + an) 



r(i + n + f)r(- 



-m — an] 



As was established in [3], there exists a set of irrational numbers £, which is uncountable, dense, 
and has Lebesgue measure zero, such that for all irrational a ^ C we have 



p{x) = x-^-^Y.ll^rn,n+iX-", if a G (0,1), (4) 

n>0 m>0 

p{x) = a:"''-i^^a„„„x'"+"", ifaG(l,2). (5) 



n>0 m>0 



(iv) Theorem 2 in [6] states that the previous result cannot be substantially improved: There exists an 
uncountable dense subset £ C £, such that for all a G £ and almost all p the series in (4), (5) do 
not converge absolutely for all x > 0. 
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This situation is clearly deficient, since we do not have a useful expression for p{x) if a G £ and we 
do not know whether the series (4), (5) converge ii a E C \ C Also, determining whether a belongs to 
£ or £ is problematic, as this would require full knowledge of the continued fraction representation of 
a (see [3, Proposition 1] and [6, Proposition 1]). Since it is impossible to have absolute convergence of 
the series for all irrational a, our only remaining possibility is to try to find some sort of "conditional" 
convergence. Absolute convergence implies that the order of summation does not matter. In the absence 
of absolute convergence, the way in which we compute the partial sums of the double series becomes 
very important. Our main result shows that the right way to compute the partial sums in (4), (5) is 
over triangles {m + an < C : m>0,n>0}c Z^, and that one can always find an increasing sequence 
Ck -> +00 (which depends crucially on the arithmetic properties of a), such that the partial sums will 
converge to p{x). 



2 Main result 



For X G M, let [x] denote the largest integer not greater than x and let {x} := x- 
part of X. The continued fraction representation (see [4]) is defined as 



\x] denote the fractional 



X 



[ao;ai,a2. 



oo + 



ai 



02 



where ao G Z and Oj G N for z > 1. The coefficients of the continued fraction can be computed recursively 
as follows: Define Xi := {x} and Xj+i := i > 1, then qq = [x] and = i > 1. For x ^ Q 

the continued fraction has infinitely many terms; truncating it after n steps results in a rational number 
Pn/ln '■= [oo; fli, 0-2, On], which provides the best rational approximation to x (among all rational 
numbers with denominators not greater than g„) and is called the n-th convergent (see [4, Theorem 17]). 
The numerators and denominators of convergents are known to satisfy the two-term recurrence relation 



anPn-l 
O'nQn—l 



■Pn-2, 
(ln-2i 



Q-l - 



0, 



P-2 

q-2 -- 



0, 



The following theorem is our main result in this note. 
Theorem 1. Assume that a ^ Q. Then for all x > 

lim b. 



p{x) 



X 



m+l+a(n+i)<ijj, 
m>0,n>0 



x'^P-^ lim 



E 



m+an 



m.+l+Q(n+i)<ijfe 
m>0, n>0 



(6) 



if a G (0,1), 



if a G (1,2), 



(7) 



where am^n one? are defined by (2) and (3) and = qkC^/a) is the denominator of the k-th convergent 
for 2/ a. 
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To obtain the proof of Theorem 1 one would need to modify only one step in the proof of [3, Theorem 
2]. Instead of reproducing the long and technical proof of the latter result in full entirety, we only explain 
here which step has to be modified. 

Assume that a G (1,2) and that a ^ Q. Combining formulas (17), (21), (24) and (26) in [3] we find 
that for every integer k>l 



m+an 



+ efc(x) 



(8) 



m+o(n+i)<fc 



m>0, n>0 

where the error term ek{x) can be bounded from above as follows 

k 



/I 

\ek{x)\ < {A{1 + x))^ X e-^^i"^'^) x JJ sec I — 



(9) 



for some constants A > and e > (which can depend on {a,p) but not on x or k). The main step 
in the proof [3, Theorem 2] is to construct a set £ C M \ Q (which is dense, uncountable, of Lebesgue 
measure zero), such that for all a ^ £ we have an upper bound 



n 

1=1 



sec 



< B3\ k>l. 



for some constant B = B{a). Combining (10) with (9) implies that efc(x) — > as 
us 



P[X) 



x^'P-^ lim 



E 



m+ctn 



(10) 

-oo, which gives 
(11) 



m+a{n+^)<k 
m>0, n>0 



To obtain the series representation (5) we use formula (2) and Proposition 1 and Lemma 1 in [3] and 
check that for a ^ L the series (5) converges absolutely, therefore the order of summation in (11) does 
not matter, and the sum can be rewritten in the form (5). 

The assumption a ^ £ is crucial for deriving inequality (10) and formula (11). The upper bound 
(10) is not true for all irrational a, and for a suitable a the product in (10) can grow arbitrarily fast. 
The following example illustrates this phenomenon. 

Example 1: We define r via its continued fraction representation 

r=[ao;al,a2,...] = [l,2,2^2los^...] 

2 

where the coefficients a„ are defined as a„+i = 2'^", n > — 1 and the numerators p„ and the denominators 
g„ of the n-th convergent are computed recursively via (6). We find the first few terms of p„ and g„ to 
be 



[P0,Pl,P2, 



[1,3,49, 



[go,gi,g2, 



[1,2,33, 



Let us take a = 2/t 1.34693878.... Since a„ are even integers for n > 1, it follows from (6) that pn 
are all odd numbers, so that we can write pn = 2r„ + 1 for some integers r„. From Theorem 13 in [4] we 
know that 



\qnT -Pn\ < 



1 



< 



1 



Q'n+l f^n+l^ln + In-l O'n+1 
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therefore |g„/a — r„ — 1/2| < 2 Using the inequahty | cos(7rx)| < tt\x — 1/2|, we conclude that 



sec 



a 



> 



cos 



(^(^--n))| -^|^-r„-i 



> 



IT 



therefore for this particular choice of a the full sequence of products in (10) cannot be bounded from 
above by any exponential function of k. 

As the previous example demonstrates, one can construct a so that the full sequence of products in 
(10) can grow arbitrarily fast. What is remarkable is that for any irrational a there is always an infinite 
subsequence of these products which are bounded from above by an exponential function of the index. 



Lemma 1. Assume that r ^ Q and r > 0. There exists a constant C = C{t) > such that for all k > 1 

(12) 



JJ |sec(7r/r)| < CG'^K 



1=1 

where qk = qki'^r) is the denominator of the k-th convergent for 2r. 

Proof. We use the following result (see [1, Lemma 4] or [7, Lemma 4]): for any f3 > 0, f3 ^ 

i 

where Qk = qk{P)- This is equivalent to 



lim — V ln(2|sin(7r//3)|) = 0, 



lim 



\} |sm(7r//3)| 



Ik 



which implies the existence of a constant C = C{/3) > 0, such that for all A; > 1 



I sm(7r/a;)| 



Using the identity sin(7r//3) = 2 sin(7r/^) cos(7r/^) we conclude 



2 ^ ^^"V" 2 



t, |cos(7r/f)| fj^ 



(13) 



Taking /3 = 2r we obtain (12). 



□ 



Combining (12) with (8) and (9) gives us the statement of Theorem 1 for a e (1,2). The proof for 
a G (0, 1) can be obtained in a similar way. 
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